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Abstract 

In the framework of the model inter-comparison study - Asia Phase II (MICS2), 
where eight models are compared over East Asia, this paper studies the influence of 
different parameterizations used in the aerosol module on the aerosol concentrations 
of sulfate and nitrate in PMio- 

An intracomparison of aerosol concentrations is done for March 2001 using dif- 
ferent configurations of the aerosol module of one of the model used for the inter- 
comparison. Single modifications of a reference setup for model configurations are 
performed and compared to a reference case. These modifications concern the size 
distribution, i.e. the number of sections, and physical processes, i.e. coagulation, 
condensation/evaporation, cloud chemistry, heterogeneous reactions and sea-salt 
emissions. 

Comparing monthly averaged concentrations at different stations, the importance 
of each parameterization is first assessed. It is found that sulfate concentrations 
are little sensitive to sea-salt emissions and to whether condensation is computed 
dynamically or by assuming thermodynamic equilibrium. Nitrate concentrations are 
little sensitive to cloud chemistry. However, a very high sensitivity to heterogeneous 
reactions is observed. 

Thereafter, the variability of the aerosol concentrations to the use of different 
chemistry transport models (CTMs) and the variability to the use of different pa- 
rameterizations in the aerosol module are compared. For sulfate, the variability to 
the use of different parameterizations in the aerosol module is lower than the vari- 
ability to the use of different CTMs. However, for nitrate, for monthly averaged 
concentrations averaged over four stations, these two variabilities have the same 
order of magnitude. 
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1 Introduction 



The model inter-c omparison study - Asia Phase II (MICS2), which follows 
the MICS Phase I (jCarmichael et all |2002| ). aims at comparing transport and 
deposition of sulfur, nitrogen compounds, ozone and aerosols in East Asia. 
Eight models (Mi, Mg) are compared for the following four periods: March 
2001, July 2001, December 2001 and March 2002. 

The eight models employed for the intercomparison use the same emission 
data, meteorological data, boundary conditions. Note that meteorological data 
that are different from the standard MICS meteorological file have sometimes 
been used, and that the size of the domain of study may differ from domain to 
domain. Apart from input data and the choice of the computational domain, 
the variability among the models is related to numerical schemes and param- 
eterizations used for transport, diffusion, deposition, scavenging and chemical 
mechanisms. Moreover, the eight models use different parameterizations in the 
aerosol module. 



This paper aims investigating to which extent the use of different parame- 
terizations in the aerosol module impact the aerosol concentrations of sulfate 
and nitrate. An intracomparison of aerosol concentrations is done in March 
2001 using different configurations of the aerosol module of one of the model 
used for the inte rcomparison. The intracomparison i s done using the model 



M«. i.e. Polair3D (IBoutahar et all (120041 ) . iFahev et all (120051 ) . iTombette et al. 



(120051 ). ISartelet and Hayamil (120041 )). that has been design e d for such multi- 



configuration works (see for instance iMallet and Sportissd (120061) ). Informa- 
tions a bout the domain of study and input data may be found in lCarmichael et al 
(l2007h . 



Hayami et al.l (120071 ) have focused on the intercomparison for aerosols, more 
especially on inorganic components: sulfate, ammonium and nitrate in PMio- 
They show that models' predictions tend to be closer to observations for sul- 
fate t han for nitrat e . Sulf ate tends to be slightly underestimated in March 
2001. lHayami et al.l (120071 ) found that the amount of sulfate is lower for the 
model M 8 than for M 3 , M 5 and M 7 , although these four models have com- 
parable amount of total sulfur. This can probably be explained by the omis- 
sion of aqueous chemistry in Mg. Although large discrepancies exist between 
models, especially for nitrate concentrations, total nitrate is consistently un- 
derestimated by most models. The nitrate concentrations predicted by Mg are 
comparable to other models, but the total nitrate (nitrate + HNO3) is slightly 
lower, which can be due to the omission of heterogeneous reactions. In Mg, 
most of the total nitrate is in the particulate phase. 
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This article is organised as follows. First, the different aerosol modules used in 
the different models that participated to the MICS study are presented with 
an accent on the aerosol module used in Polair3D. Then, the sensitivity to the 
aerosol module is studied and compared to the sensitivity to the chemistry 
transport model for monthly averaged concentrations at EANET (Acid De- 
position Monitoring Network in East Asia) stations and at Fukue, a remote 
site between Japan and China. 



2 The aerosol modules 



2. 1 The different physical processes 



In the aerosol module of eulerian models, the si z e dist ri bution is o f ten m odeled 
using a sectional ap proach (e.g. IPilinis et al.l ( 2000[). Ijacobson or a 



modal approach (e.g. lWhitby and McMurryl (119971 ). ISartelet et al.l (120061 )). In 
the sectional approach, the number of sections varies from study to study, 
usually from 2 to 15. In the modal approach, typically from 2 to 4 modes are 
used. 



The aerosol composition and distribution are influenced by different phys- 
ical processes such as condensation/ evaporation, which is thought as one 
of the major processes influencing aerosol composition, heterogeneous reac- 
tions on the surface of particles, coagulation, nucleation, cloud chemistry. 
Three approaches may be used to model condensat ion/ evaporation processes 
(ICapaldo et ail fj200oh . iDebrv and Sportissel (l2006t )) 



a dynamic approach, i.e. the mass transfer between gas and aerosol phases 
is explicitly taken into account, 

a full-equilibrium approach, i.e. the dynamic modeling is replaced with 
an assumption of thermodynamic equilibrium between the gas and aerosol 
phases, 

an hybrid approach, i.e. full equilibrium is assumed for fine aerosols while 
the dynamic approach is used for coarse aerosols. 



2.2 The aerosol modules in MICS models 



The main differences of the different aerosol modules used in the MICS models 
can be classified as differences in the size distribution, i.e. modal versus size re- 
solved, the number of sections or modes, and differences in physical processes, 
i.e. coagulation, condensation/evaporation, cloud chemistry, heterogeneous re- 
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actions. Although input data, such as boundary conditions and emissions, are 
common to all models in the MICS study, sea-salt emissions are not given as 
input data but they are parameterized in some models. These differences are 
summarized in Table [TJ 



2.3 PolairSD 



he se nsitivity study is carried out with the model M 8 , i.e. Polair3D (B outahar et al 



(2004)). The aerosol module of Polair3D is fully described in iDebry et al. 



(120061 ). Aerosols are composed of black carbon, dust, fiv e inorganics ( sulfat e, 
nitrate, ammonium, sodium, chloride) and eight organics ( Schell et al.l ( 2001 )). 
The number of sections used in the modeling is fixed by the user. The sections' 
diameters are log-distributed in the diameter size-range considered, 0.01/mi 
to 10/im in this study. The bounds of the sections are 0.01, 0.02, 0.04, 0.08, 
0.16, 0.32, 0.63 , 1.26, 2.51, 5.01, 10 . 00//m . Brownian coagulation is modeled 
as described in Debrv and Sportissel (2007 ). and condensation/evaporation as 
described in Debry and Sportissel ( 20061 ) with a moving sectional scheme. The 
three approaches, the full equilibrium approach, the dynamic approach or the 
hybrid approach, m ay be used. For cloud che mistry, the variable size resolved 
module (VSRM) of iFahey and Pandid (120031 ) can be switched on or of f. Het- 
erogen eous reactions on the surface of particles are modeled following [Jacob 
(l2000h 

H0 2 — >0-5H 2 O 2 (Rl) 



NO, 



0-5HONO + 0-5HNO, 



(R2) 



NO, 



HNO, 



(R3) 



N 2 5 



2HNO, 



(R4) 



They are considered irreversible processes with a first order reaction rate. 



The parameterizations of sea-salt emissions and vertical diffusion a re prepro- 

cessed. Sea-salt emissions are computed with the paramet erization of iMonahan et al 
( 119861 ) generalized to varying relative humidity following IZhang et al.l (120 05'). 
Vertic al diffusion is computed with the parameterization of iTroen and Mahrt 
jl986h . 



2.4 The sensitivity study 



In this sensitivity study, only the concentrations of sulfate and nitrate are 
studied. The sensitivity tests concern the size distribution, i.e. the number of 
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No 
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No 
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Table 1 
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sections, and physical processes, i.e. coagulation, condensation/evaporation, 
cloud chemistry, heterogeneous reactions and sea-salt emissions. Sensitivity to 
nucleation is not tested because nucleation would involve particles of size be- 
tween 0.001//m to 0.01/xm, a size range smaller than the one considered in this 
study (0.01/xm to 10/j.m). The sensitivity to whether condensation/evaporation 
is treated dynamically or assuming thermodynamic equilibrium is also studied. 
However, the sensitivity to the ther modynamic model i s not tested, although 
it may not be negligi ble as s hown in IZhang et al.l (119991 ) . Heterogeneous reac- 
tions as described in I Jacob! (120001 ) are coded in Polair3D. Other models used 
in MICS may however take into account different heterogeneous reactions, or 
similar ones but with a different reaction rate. 



In the default run Ri, ten sections are used, condensation/evaporation pro- 
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No 
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Table 2 

Summary of the different options used in the runs performed for the sensitivity 
study (these configurations refer to model M$). 

cesses are modeled with the full equilibrium approach, coagulation is taken 
into account, cloud chemistry, sea-salt emissions and heterogeneous reactions 
are ignored. In each of the runs R2 to R7, one option differs from the run 
R\. For example, in the run i? 2 , only three sections are used. A summary of 
the different runs is presented in Table |2j In the run i? 3 , the hybrid approach 
is used with a cutoff diameter of 0.6/im: the finest 6 sections are computed 
with the full equilibrium approach, while the dynamic approach is used for 
the coarsest 4 sections. To show that the aerosol concentrations are not only 
influenced by the parameterizations in the aerosol module, a run is done with 
a different parametri z ation of the vertical diffusion. The sensitivity study of 



Mallet and Sportissd (120061 ) shows a great sensitivity to the parametrization 
of the vertical distribution for ozone. Therefore, the parameterization of Louis 
is used instead of the parametrization of Troen-Mahrt in the run R§. 

The sensitivity of sulfate and nitrate concentrations to the different options is 
now assessed by comparing monthly averaged concentrations. 



3 Monthly averaged concentrations 



Monthly averaged concentrations of sulfate and nitrate were measured at some 
of the EANET (Acid Deposition Monitoring Network in East Asia) stations 
( IHayami et all 120071 ) . The sensitivity to the aerosol module is assessed by com- 
paring the different runs at 23 EANET stations and at a remote site betwee n 
Japan and China, Fukue, a station operated by CRIEPI (IHayami et all 120071 ) . 
The location of the stations is given in Figure [TJ Measurements for aerosols, 
that are in the computational domain of the eight models, are available at only 
three of the EANET stations (Terelj -Mongolia, Kanghwa and Imsil -Republic 
of Korea) and at Fukue. 
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Hedo 



Fig. 1. Location of the stations used in the sensitivity study 
3.1 Sensitivity to the aerosol module 



To quantify the sensitivity of the different aerosol options, the normalized 
mean absolute error (NMAE) between Ri and Ri (i = 2, ...,8) and the nor- 
malized mean b ias (NMB) are computed for sulfate and nitrate, as follows 
(lYu et all fl2006h ) 



NMAE = ^ fc=1 7' fc lM . 100% 

NMB = Efc=1 ffi'* ~ Cl ' k) ■ 100% (1) 

where N is the number of stations (23 EANET stations + Fukue), C^k (i = 
1, 8, k — 1, N) represents the concentration of sulfate or nitrate at the 
station k for the run For each run Ri, the higher the NMAE and NMB are, 
the more sensitive the results are to the option i. A summary can be found in 
Table El 



As shown in Table El sulfate concentrations are little sensitive to options 3 
and 7, ie. to the hybrid option and to sea-salt emissions, with a NMAE under 
2%. Sulfate is little sensitive to sea-salt emissions, because sea-salt sulfate is 
not considered in the sea-salt emissions here. Sea-salt emissions are assumed 
to be made exclusively of chloride and sodium. 

Higher sensitivity of sulfate is observed for options 2 and 4, i.e. for the size 
distribution and for coagulation with a NMAE of 6% and 11% respectively. 
Sulfate concentrations are equally sensitive to options 5 and 6, i.e. to het- 
erogeneous reactions and cloud chemistry, with a NMAE around 15%-17%. 
Although the formation of sulfate is often thought to be dominated by aque- 
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R3 
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-10% 
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-15% 


R5 
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-14% 
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99% 


Re 


17% 


17% 


2% 


-1% 


R 7 


3% 


3% 


10% 


10% 


Rg 


34% 


34% 


71% 


71% 



Table 3 

NMAE and NMB between R x and R4 (i = 2, 
trations at 23 EANET stations and at Fukue. 



, 8) , for monthly averaged concen- 



PRECIP(mm/month) Z: srfc. MAR2001 




50 100 150 200 250 500 550 



Fig. 2. Precipitations in mm/month for March 2001. 



ous production (e.g. Seinfeld and Pandisl ( 1998 )). the sensitivity of sulfate to 
cloud chemistry is only 17% here. However, this sensitivity varies strongly 
with the location of stations. Figure [2] shows the total amount of precipita- 
tions for March 2001. Precipitations are low over Mongolia and central China, 
but high over south China and Japan. Accordingly, the sensitivity of sulfate 
to cloud chemistry is only 1% at Terelj in Mongolia (see Figured]), whereas it 
is as high as 71% and 62% in Xiaoping and Hedo. Note that higher sensitivity 
to cloud chemistry would have been observed in July, when precipations are 
higher than in March. 



NMAE tends to be higher for nitrate than for sulfate concentrations, except for 
the option 6, i.e. for cloud chemistry. Nitrate is more sensitive to the options of 
the aerosol module than sulfate, explaining the larger discrepancies observed 
between the models M\, ...,M 8 for nitrate than for sulfate. The effects of the 
size distribution, the hybrid scheme and sea-salt emissions are stronger on 
nitrate than sulfate because they alter thermodynamic equilibrium. Sulfate is 
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less affected because of its low volatility. 

For nitrate, the highest sensitivity is observed for option 5, i.e. heterogeneous 
reactions, with a NMAE as high as 99%. Only cloud chemistry (option 6) 
does not significantly modify nitrate concentrations. The effects of the size 
distribution, the hybrid scheme, coagulation and sea-salt emissions (options 
2, 3, 4 and 7) are of similar strength in terms of NMAE with a NMAE around 
10%-15%. 

The positive NMB of the run Rq for sulfate illustrates that the concentration 
of sulfate increases due to cloud chemistry as expected by the oxidation of 
dissolved S0 2 into sulfate. 

Nitrate concentrations increase by heterogeneous reactions as shown by the 
high positive NMB of 99% for the run R 5 , whereas sulfate concentrations de- 
crease although the NMB —14% is not as high in absolute value as for nitrate. 
As deduced from the set of reactions IR 1\ IR21 IR3I and IR4[ heterogeneous 
reactions lead to higher HNO3 and H2O2 concentrations and to lower HO2 
concentrations. Higher H2O2 concentrations in R§ would increase sulfate con- 
centrations in the real atmosphere. However, this does not occur in the model 
runs because R 5 excludes aqueous chemistry. A more detailed comparison of 
the runs Ri and R^ shows that monthly averaged OH and O3 concentrations 
are also lower when heterogeneous reactions are taken into account. Because 
of lower OH concentrations, oxidation of SO2 by OH is lower resulting in 
lower sulfate concentrations. Lower sulfate may in turn lead to higher nitrate, 
which is then required to neutralize ammonia. Higher nitrate concentrations 
may also be a consequence of the higher HNO3 concentrations that condense 
on particles to form ammonium-nitrate. 

Coagulation is more efficient for small particles, limiting their concentrations. 
Figure [3] shows the size distribution of sulfate and nitrate monthly averaged 
concentrations, averaged over all stations. When coagulation is not taken into 
account, the distribution is centered at smaller diameters than when coagula- 
tion is taken into account. 

For both sulfate and nitrate, a high sensitivity to the parameterization of the 
vertical diffusion (run R§) is observed. For sulfate, this sensitivity is much 
higher than the sensitivity to any of the options tested in the aerosol module. 
For nitrate, this sensitivity is very high as well, although not as high as the 
sensitivity to heterogeneous reactions. 
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Fig. 3. Size distribution of sulfate and 
averaged over all stations. 

3.2 Comparison to the sensitivity to 



Rl (default) 
R4 (no coag) 




Diameter (|im) 

nitrate monthly averaged concentrations 
the chemistry transport model 



To compare the sensitivity to the aerosol module to the sensitivity to the 
chemistry transport model, the NMAE is computed for stations where mea- 
surements are available (Terelj, Kanghwa, Imsil, Fukue) as follows 



NMAE = fc=ll A 7" — • 100% (2) 



where M = 4 is the number of stations where measurements are available, Ok 
(k = 1, ...,N) is the concentration of sulfate or nitrate at station k. Table H] 
shows the mean NMAE for runs R\ to Rj as well as the smallest (Min) and 
largest (Max) value of NMAE obtained for these runs. These values Min, Max 
and mean of NMAE are computed not only for runs Ri to R7 but also for 
models M\ to Ms, in order to compare the sensitivity of the aerosol module 
(runs Ri to R7) to the sensitivity of the chemistry transport model (models 
Mi to M 8 ). For sulfate, all eight models are included when computing the 
mean NMAE. For nitrate, the model M3 is not included because it did not 
provide results. Furthermore, because the NMAE observed with the model M4 
is one order of magnitude higher than other models for nitrate, results with 
and without taking M 4 into account are presented. 

Values of NMAE are higher for nitrate than sulfate, confirming again the 
higher sensitivity of nitrate. 

For sulfate, variability is larger for models Mi to M 8 than for runs Ri to 
R7. The NMAE varies by only 2% for runs Ri to R7 whereas it varies by as 
much as 8% for models Mi to M 8 . For nitrate, the NMAE varies by 65% for 
runs Ri to R-j whereas it varies by as much as 175% for models Mi to M 8 . 
However, without taking into account the model M4 for nitrate, the NMAE 
varies by 56% for models Mi to Ms, that is the variations are of the same 
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Sulfate 



Nitrate 





Min 


Mean 


Max Min 


Mean 


Max 


Ri — R-j 


12% 


13% 


14% 20% 


34% 


85% 


Mi - M 8 


10% 


14% 


18% 13% 


57% 


188% 


Mi - M 8 (no M 4 ) 






13% 


35% 


69% 



Table 4 

Min, mean and Max of NMAE associated to runs Ri to Rj and to models Mi to 
Mg for monthly averaged concentrations. 

order of magnitude as those of runs Ri to R-?. In other words, for nitrate, the 
variability within the aerosol module (Ri to R7) is as large as the variability 
between models (Mi to M 8 ). 

Apart from the uncertainty in aerosol parameterizations, the variability among 
models Mi to Mg is linked to numerical schemes and parameterizations used 
for transport, diffusion, deposition, scavenging, chemical mechanism, choice 
of computational domain and input data. Although input data are the same 
for all models for emission and boundary conditions, meteorological data that 
are different from the standard MICS meteorological file have sometimes been 
used. For example, a difference in t he relative humidity field used by models 



as high as 20% has been observed (jHozumi et all 120051 ). Such variability in 



relative humidity contributes to variability in aerosol concentrations. 

Variability linked to numerical schemes and parameterizations outside the 
aerosol module is important especially for sulfate for which the variability 
to parameterizations in the aerosol module is not so high. For example, the 
sensitivity to the parameterization of the vertical diffusion has shown to be 
higher than the sensitivity to any options of the aerosol module for sulfate. 

Because only uncertainties in aerosol parameterizations are considered for the 
runs Ri to R-j, uncertainties observed for models Mi to M 8 should be larger. 
This is true for sulfate. However, for nitrate, if the results of the model M 4 are 
not considered, the NMAE varies as much for models M\ to M 8 as for runs 
R\ to Rj. This strong variability of nitrate to the aerosol module stresses the 
difficulties in modeling nitrate concentrations accurately. 



4 Conclusion 



For March 2001, the sensitivity to the aerosol module of the model Mg, Po- 
lair3D, is assessed and compared to the sensitivity to the chemistry trans- 
port model. Because in the MICS comparison, all models are assumed to 
use the same input data, the variability between the different models is due 
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to differences in physical parameterizations, differences in numerical schemes 
and differences in the chemical mechanism. Concerning the sensitivity to the 
aerosol module, only the sensitivity to the size distribution (number of sec- 
tions) and the sensitivity to physical parameterizations (coagulation, whether 
condensation is solved assuming full equilibrium or using an hybrid scheme, 
cloud chemistry, heterogeneous reactions, and sea-salt emissions) are consid- 
ered. However, there is also a sensitivity to numerical al gorithms for s imula tion 
of condensation/evaporation for example as shown by Zhang et al.l ( 2004 ). 



To assess the sensitivity to the aerosol module and to the chemistry trans- 
port model, monthly averaged concentrations over 24 stations are computed. 
Sulfate concentrations are little sensitive to the hybrid scheme and to sea-salt 
emissions, with a normalized mean absolute error (NMAE) and a normalized 
mean bias under 2%. Nitrate concentrations are little sensitive to cloud chem- 
istry. For sulfate, the sensitivity to cloud chemistry is important although not 
dominant with a NMAE of about 15%. For sulfate and especially for nitrate, 
high sensitivity to heterogeneous reactions is observed. For nitrate the sensi- 
tivity to heterogeneous reactions is particularly high, with a NMAE as high 
as 99% for monthly averaged concentrations. 



To compare the sensitivity to the chemistry transport model to the sensi- 
tivity to the aerosol module, the minimum, the mean and the maximum of 
the NMAE are computed for R iy i.e. for runs associated to sensitivity tests 
on the aerosol module, and for Mj, i.e. for models participating in the MICS 
comparison. The values of NMAE are higher for nitrate than sulfate, con- 
firming the higher sensitivity of nitrate. For monthly averaged concentrations 
at all EANET stations for nitrate, the variations of the NMAE for Ri are of 
the same order as the variations for Mj, suggesting a very high sensitivity to 
the aerosol module. For sulfate, the variations of the NMAE for Ri are much 
smaller than the variations of the NMAE for Mj (2% against 8%), suggesting 
that the sensitivity to the aerosol module is less strong than the sensitivity 
to the chemistry transport model. However, if R s , i.e. the sensitivity to the 
parameterization of the vertical diffusion, is added in the computation of the 
variations of the NMAE for Ri, then the variations of the NMAE for Ri are 
almost of the same order as the variations for M, (6% against 8%). 



The high variability observed for nitrate stresses the difficulties in modeling 
nitrate concentrations accurately. Because nitrate concentrations are highly 
sensitive to heterogeneous reactions, better knowledge of reaction probabil- 
ities for example, some of which may be expressed as a function of a eroso l 
composition, temperature and relative humidity (e.g. lEvans and Jacob! (120051 ) 
for IR 4D . may help to improve the models accuracy. 
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